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We study the long-time behavior of underdamped Brownian particle moving through a viscous 
medium and in a systematic potential, when it is subjected to a space-dependent high-frequency peri- 
odic force. When the frequency is very large, much larger than all other relevant system-frequencies, 
there is a Kapitsa time- window wherein the effect of frequency dependent forcing can be replaced by 
a static effective potential. Our new analysis includes the case when the forcing, in addition to being 
frequency-dependent, is space-dependent as well. The results of the Kapitsa analysis then lead to 
additional contributions to the effective potential. These are applied to the numerical calculation of 
the diffusion coefficient (D) for a Brownian particle moving in a periodic potential. Presented are 
numerical results, which are in excellent agreement with theoretical predictions and which indicate 
a significant enhancement of D due to the space-dependent forcing terms. In addition we study the 
transport property (current) of underdamped Brownian particles in a ratchet potential. 
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I. INTRODUCTION 

Brownian motion in periodic structures has various 
applications to condensed matter physics, nanotechnol- 
ogy and molecular biology Q, |^ 1^. Adding noise to 
deterministic nonlinear dynamics has led to interesting 
and important phenomena such as stochastic resonance 
0, Brownian motors and chaotic ratchet transport 0, 
resonant activation j^, noise induced phase transition 
0, 111, etc.. Thermal diffusion of a Brownian particle 
which we will discuss here is of great interest in numerous 
other contexts, namely Josephson junction 0, rotating 
dipoles in external fields ^3 , superionic conductors ■, 
synchronization phenomena |l2l| . diffusion on crystal 
surfaces particle separation by electrophoresis [T^ . 
and biophysical processes such as intracellular transport 

El. 

In this paper we focus on the underdamped motion of a 
Brownian particle which feels viscous forces and random 
impulses from the surrounding medium and is confined 
by a potential well. Our primary interest is to study 
the effect of an externally applied position dependent 
driving force that is periodic in time. The frequency is 
much larger than all other relevant frequencies of the 
system. Hence we can apply the usual Kapitsa analysis 
for high frequency oscillating fields It has earlier 

been shown that on time scales larger than the period of 
perturbation the dynamics is equivalent to one in which 
the periodic perturbation can be replaced by a time 
independent effective potential [T^ . ITq . Ref . 01 treats 
the overdamped Brownian motion, whereas we deal with 
underdamped motion of the Brownian particle. Dutta et 
al. E3 have based their analysis on the Fokker-Planck 
equation approach. Here we provide an alternative 
derivation of the main results through the Langevin 
dynamics, which is more straightforward. 
Further we extend Kapitsa's analysis by using the addi- 



tional contributions to the effective potential arising from 
the space dependent periodic force for the calculation 
of the thermal diffusion coefficient 0, |23| . The central 
result of our present paper is that the effective diffusion 
coefficient of an underdamped Brownian particle in a 
periodic potential in presence of an externally applied 
space-dependent periodic force can become arbitrarily 
larger than in the presence of thermal noise alone. 
Furthermore, certain features of the ratchet mechanism 
|2ll | are relevant in this context in terms of transport 
properties (currents) . 

With the preceding background the paper is organized as 
follows. In Sec. II we introduce the model, the numerical 
scheme and the basic quantities of interest, namely, the 
effective potential, the effective diffusion coefficient, and 
the average particle current. In Sec. Ill we develop the 
necessary formalism to address rapidly periodic drive 
and arrive at the perturbative effective potential. In 
Sec. IV we discuss the numerical results on the effective 
diffusion coefficient and average particle current. The 
summary remarks, discussion and conclusion of our 
findings are presented in Sec.V. 

II. MODEL 

The stated Brownian dynamics is governed by the 
Langevin equation 

mi = -■^x- —U{x)+F{x,t)+7]{t), (1) 

where m is the mass of the Brownian particle, 7 is the 
friction coefficient, U(x) is the confining potential, F(x,t) 
is the periodic driving force with a period r and F(x,t) = 
F(x,t-|-T). Thermal fiuctuations are modeled by the zero 
mean 5-correlated white noise i]{t) i.e. < ri{t) > = and 
< r]{t)Tj{t') >= 2-fp-^6{t - t'), where /3 = [keT)-^, kg 
being the Boltzmann constant and T is the temperature. 
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Our goal is to show that, on time scales larger than r, 
the dynamics can be mapped onto a modified Langevin 
dynamics in which the periodic forcing is absent but the 
potential U(x) can be replaced by a suitable effective po- 
tential Ueff{x). The methdology we follow is based on 
Kapitsa's treatment for high frequency oscillating field 
in parametric oscillations |22 | . We derive the form of the 
effective potential upto second order in ^ (expansion pa- 
rameter which is basically inverse of the square of the 
oscillating frequency) in Sec. III. 

We also study the transport properties and diffusion 
coefficient for ratchet like systems. The corresponding 
Langevin dynamics is governed by the equation 

-V'ix) + A{x) cos(m) 4- V27^^(0, (2) 

where V(x) is the periodic potential with period L i.e. 
V(x) — V(x-I-L) and the prime denotes 1st derivative of 
V(x) with respect to x. In our case V{x) = — sin(a;) — 
fism{2x) and ji = j throughout this work. We define 



1.5 




-1.5 I ' ' ' 1 

-10 -5 5 10 

X 

FIG. 1: The periodic potential V(x). 

the above dynamics as the original dynamics. Following 
Kapitsa's treatment we derive in the sequel an effective 
potential for which the dynamics is governed by the equa- 
tion 

mi + -fx = ~V^ff{x) + y/2-fkBTi]{t). (3) 

The first basic quantity of interest in the transport 
process is the average particle current defined as 

< X >= lim — — . (4) 

t—toc t 

The other quantity of important interest is the effective 
diffusion coefficient 

i^=lim<"'^^)^:<"^^) (5) 

Exact analytical results for D are known for two special 
cases. First in the absence of the periodic potential we 
have the famous Einstein's relation D = Second is 

7 

the case in which the external field A{x) cos(fit) is ab- 
sent. Calculation of the diffusion coefficient in the pres- 
ence of both externally applied space-dependent periodic 
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FIG. 2: The effective potential 14// (x). 

force and arbitrary periodic potential is not analytically 
possible. Hence we solve Eq.(2) and Eq.(3) numerically. 
To compute x(t) we need to transform Eq.(2) and Eq.(3) 
into dimensionless Langevin equations as follows: 

x + bx ^ -V'{x) + a{x) cos(wi) + y/2bDofj (6) 
x + bx = -Kffi^) + V'^bDoV, (7) 

where Eq.(6) denotes original dynamics, whereas Eq.(7) 
is for effective dynamics. In writing the above equations 
we define the following dimensionless quantities: 

1. the rescaled friction coefficient b = -^tq which is 
the ratio of the two characteristic time scales, rg 
and Ty. To is determined from Newton's equation 
mx — ~V'(x), by balancing the two forces = 

^ yielding = whereas ry = 

2. the rescaled potential V{x) = -^^j assuming pe- 
riod 1 and barrier height AV = 1, 

3. the rescaled force strength a ~ with the fre- 
quency uj = r^To, 

4. the rescaled white noise forces fi{t) obeying 

< f][t)f}{t') >= 5{t—t') with rescaled noise intensity 

III. EFFECTIVE POTENTIAL 

As mentioned earlier our focus is on the result that 
on time scales larger than the period of perturbation, 
the dynamics is equivalent to one in which the time 
dependent periodic perturbation can be rep laced by a 
time independent effective potential l23| . In Sarkar 
and Dattagupta it has been shown in great detail 
that the expression of the effective potential does not 
alter in the presence of noise. We presume and verify 
that this result is true even when the forcing term 
is space dependent. Further, unlike we treat the 
underdamped case from which all the results of can 
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be obtained as a limit. 

A. First Order Correction 

It is evident from the nature of the field in which the 
particle moves that it will traverse a smooth path and at 
the same time will execute small noisy fluctuations about 
that path. Accordingly, wc represent the function x(t) as 



x{t)=X{t)+at), 



where X(t) is a slow variable and ^(X, t) is a fast variable. 
The following transformations then follow: 



i = X + ^{X,t) 
X = X+'i{X,t) 



d_ 

dx 



_l d_ 



(9) 



Now setting the noise term to zero and putting the above 
transformations in Eq.(l) we obtain 



m 



(x{t)+'i{X,t)) = -f(^X{t)+iiX,t))+F{X + tt) 

d 



1 



1 + CdX 



u{x + o. 



(10) 



To find the effective potential experienced by the particle 
correct to first order, we perform a Taylor series expan- 
sion of Eq.(lO) upto first order. Thus, 

mX{t)+mi{X,t) = --fX{t) --fi{X,t) 

^ ^ {UiX)+^U'iX)) 



l + S,dX 
+ F{X,t)+^F'iX,t). (11) 

The above Eq.(ll) involves both 'fluctuating' and 
'smooth' terms on the left and right sides which must 
be separately equal. For the fluctuating terms wc can 
simply put 



mat)+jat) = Fix,t), 



(12) 



where wc take F{X, t) = f{X) cos{ujt)+g{X) sin(wt). 
Solving Eq.(12) wc obtain, 



jnii^i + X) 



!{X) + ^g{X) coaiut) 



g{X)-^f{X))sm{ujt) 



(13) 



Since oj is large, ~ (1^0 is effective. Next combin- 
ing Eq.(12) with Eq.(ll) and then averaging over a time 
period wc finally obtain 



mX{t)+jX{t) 



dUjX) 
dX 



+ <^F'{X,t)>, (14) 



where 
< ^F'{X,t) > = -- 



4m ( uj'^ + 



— {f{X)+g\X)) 



+ ^{nX)g{X)~g'{X)f{X))\.{lh) 

With the help of Eq.(15) we can rewrite Eq.(14) as follows 

dU{X) 

d 



(8) mX{t)+-iX{t) 



dX 



1 



4m 



dX 



{r{X)+g'iX)) 



+ ^{f'iX)g{X)-g'{X)f{X)) ,(16) 



mX{t) +"/X(t) 



dUeffjX) 
dX '■ 



(17) 



with UeffiX) = U{X) + U^{X), where 



U\X) = 



1 



4m { ui'^ + 



±^{f{X)+g\X)) 



+ 
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moj 



dy{f'{y)g{y)-9'{y)f{y)) ■ (18) 



It is clear that U^{X) vanishes for spacc-indcpcndcnt 
forcing. 

B. Second Order Correction 

To find the second order correction term in the effective 
potential the transformation equations given in Eq.(9) 
have to be modified as 

X = x+ax,t)+ax,t)C{x,t) 
X = x+i+ie+it' 
i = x+e+ir+ee' + 2cc' 

d 1 d 



dx 



1 + e' + e'^ + ec" dx 



(19) 



Putting the above transformation in Eq.(l) and retaining 
terms upto seond order in ^ (0(^^)) we derive 



m( X + e + + ee' + 2ee'j = -7(^ + e + ee' + 

+U"{X)^ + [/'(X)e' + U"{X)ii' + U\X)i'^ 

+u'{x)ii" + ]f"{x)e + c/"(^)ee'] + F{x, t) 

HF'{X, t) + ^S.'F'{X, t) + \eF"{X, t). 
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Next, after performing time averaging we ultimately ob- 
tain 

mX + -iX = -U'{X) + U'{X) <£,'^ > 



1 d 
2dX 
dUeff 
dX 



{U"{X)<(' >)+<^F\X,t) > 



(20) 



with UeffiX) = U{X) + U'^{X) + U^{X), where U(X) is 
the systematic periodic potential, U^{X) is the first order 
correction to the effective potential given by Eq.(18) and 
the second order correction to the effective potential is 
given by 



U^{X) 



1 



dy(r{y)+g'\y) 



{f{X)+g\X))U"{X) 



(21) 



These results are identical to those obtained by Dutta 
and Barma |l7l | who have however employed a Fokker- 
Planck equation approach. In the next section we 
further extend this analysis by using these results for 
calculating the effective diffusion enhancement and 
transport current. 

IV. NUMERICAL SCHEME AND RESULTS 

We numerically solve the following two dimensionless 
equations of motion 



x + bx + V'{x) 
X + bx 



a{x) coa{Lut) + ^/2bDov{t) (22) 
-Uiffix) + V^bD'ovit), (23) 



with the aid of the Heun scheme which is basically the 
Runge-Kutta algorithm. Our main interest as empha- 
sized earlier is to calculate the effective diffusion coeffi- 
cient which is defined as 



= lim 

t — >oc 



« x'^{t) » - « x{t) >>^ 
2t 



(24) 



where the two brackets respectively denote averages over 
the initial conditions of position and velocity and over all 
realizations of thermal noise. We calculate this quantity 
for the orginal dynamics and effective dynamics for the 
two special cases: (a) space dependent external periodic 
force and (b) constant amplitude external periodic 
force. The results are shown in Fig. 3. In both the 
cases the enhancement of effective diffusion coefficient 
as a function of the amplitude 'a' is clearly noticeable. 
The important difference is that the enhancement is 
much more pronounced for the space dependent external 
periodic force than the constant amplitude external 
periodic force. 

It is known that there are two states of a driven 
Brownian dynamical system: the locked states, in which 
the particle stays inside one potential well, and the 




FIG. 3: The effective diffusion coefficient for original and ef- 
fective dynamics for two cases (a) space dependent (b) space 
independent external periodic force. 



running states, for which the particle runs over the 
potential barriers. The first regime is characteristic of 
a small driving force strength. When the amplitude of 
the external field is made large, running states appear 
where we can see both diffusive and regular behavior 
of the particle. We have shown above that the effective 
diffusion coefficient increases rapidly when the externally 
applied force exceeds a certain critical strength. This 
leads to the existence of an optimal "a" for the enhance- 
ment of diffusion rate. This phenomenon is reminiscent 
of Stochastic Resonance (SR) H [H [H HHH IH ■ 
So we can hereby employ the acronym "SR" to imply 
accelaration of diffusion. By this we mean that a new 
diffusion mechanism, with combined action of noise, 
spacially periodic potential and time-periodic modula- 
tion can be more effective than that of a free Brownian 
motion, since D^ff is shown to exceed unity in a large 
region around some optimal parameter regions. In these 
regions the optimal matching of the periodic force and 
the noise drive the particles up to the potential hills 
from where they scatter at the potential barriers and 
finally diffuse very quickly into wide regions. For small 
'a' it is very difficult to push them up to the potential 
hills. Therefore, in order to get the above mentioned 
diffusion enhancement we need the optimal collective 
actions of three forces: spatially periodic gradients, 
time-periodic modulation and stochastic stimulation. 
Further we should emphasize that the extra terms in 
effective potentials due to space dependence of the 
external force do indeed aid this diffusion enhancement 
mechanism. 

We have also studied the current 'J' which is defined 
as the time average of the average velocity over an 
ensemble of initial conditions. Thus it involves two 
different averages, the first is over M initial conditions, 
which we take randomly centered around the origin and 
with an initial velocity equal to zero. For fixed time tj we 
calculate the average velocity Vj = -jj X^t^i ^i(^j)- The 
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second average is over time and yields J = ^ ^j=i ^J' 
All quantities of interest are averaged over 250 different 
trajectories and 10^ periods. There are three dimen- 
sionless parameters a, b and uj defined earlier in terms 
of physical quantities. We vary the parameter 'a' and 
fix 6 = 0.1 and uj = 0.67 throughout. In Fig. 4 we have 
shown the behavior of the transport currents in case 
of space dependent external force. Initially the current 
is zero, following which it increases and peaks at some 
optimal values of 'a', then decreases with the increase of 
'a'. The following explanation will help to understand 
the behavior of current. At very low force strength, 
escape jumps between the neighbouring wells are very 
rare i.e. the average directed current is very small. The 
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FIG. 4: The current J of the inertial rocked Brownian motor 
when external force is space dependent. Inset shows the graph 
for very low strength force ('a'). 

input energy is mostly expanded into the kinetic energy 
of the intrawell motion and eventually dissipates. As 
'a' is increased further, the Brownian motor mechanism 
starts to work and some part of energy contributes to the 
net motion of the particle. Now due to inertia, the mean 
velocity increases and reaches a maximum. Then 'a' 
reaches a second threshold value above which the current 
starts to decrease because of the debilitating effect of the 
ratchet potential. The occurrence of multiple reversals 
of the directed current, as is shown in the inset of 
Fig. 4, is an interestin g fe ature of the inertial Brownian 
motor system |^ Is^. l33l Is^ Issj . The phenomenon of 
current reversals can be described by different stability 
properties of the perturbed rotating orbits of the system 

By comparing the figures (4) and (5) we can surmise 
that the current is much more substantial for the space 
dependent external force case. Extra terms in the 
effective potential arising from space dependence of the 
external force does help in increasing the current. 



V. SUMMARY AND CONCLUSIONS 

In this section we present an overview of the principal 
results of this paper. We have addressed the problem of 
underdampcd Brownian particle in a position dependent 
periodic driving force in the high frequency regime. We 
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FIG. 5; The current J of the inertial rocked Brownian motor 
when external force is space independent. Inset shows the 
graph for very low strength force ('a'). 



have then calculated the effective potential upto second 
order in the expansion parameter ^ and used these 
results to calculate the effective diffusion coefficient 
and transport current. In the high frequency regime 
the particle makes small but rapid excursions around a 
smooth path along which the motion is relatively slow. 
A systematic perturbative treatment in powers of the 
excursion amplitude shows that the first order correction 
in the effective potential exists only if the externally 
applied rapidly oscillating field is space dependent. 
This first order correction term (Eq.(18)) is the average 
kinetic energy which contributes to the work done 
against damping. The second order correction to the 
effective potential shows that a nontrivial contribution 
arises even for position independent driving. 
We have employed our derived results for the calcula- 
tion of the effective diffusion coefficient and transport 
current. We obtained the effective diffusion coefficient 
by solving both the orginal dynamics and effective dy- 
namics. We noted a giant enhancement of diffusion and 
the results arising from original and effective dynamics 
agree very well. This validates our method of calculation 
in the high frequency regime. The enhancement of 
diffusion is a result of the optimal collective actions of 
spatially periodic gradients, time periodic modulation 
and thermal noise. The enhancement is much more 
pronounced for the space dependent periodic external 
force which can be understood in terms of the extra 
terms arising in the effective potential from the space 
dependence of external force. We have analyzed the 
transport properties and the behavior of current for the 
Brownian motor mechanism and compared the currents 
for two cases: space dependent and space independent 
external forces. The current is much larger for the space 
dependent case. 

Finally we would like to emphasize once again the practi- 
cal implications of this work. The parameters specifying 
the particles enter the effective potential, which can be 
used to separate different species of Brownian particles 
by identifying the minima of the effective potential. One 
can control the diffusion rates by varying periodic spatial 
gradients and the space dependent external field. In ad- 
dition to myriad applications mentioned earlier, systems 
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described by Eq.(l) are realized for charged particles 
moving on thermal surfaces under periodic potentials 
subjected to time varying fields. Recent motivation to 
study these systems has been inspired by the theoretical 
modeling of the molecules called kinesin and myosin, 
which posess the ability to move unidirectionally along 
structural filaments of microtutubulin and actin. 
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